The landscape of the immunoglobulin repertoire in endemic pemphigus foliaceus.

Introduction: Primarily driven by autoreactive B cells, pemphigus foliaceus (PF) is an uncommon autoimmune blistering skin disease of sporadic occurrence worldwide. However, PF reaches a prevalence of 3% in the endemic areas of Brazil, the highest ever registered for any autoimmune disease, which indicates environmental factors influencing the immune response in susceptible individuals. We aimed to provide insights into the immune repertoire of patients with PF living in the endemic region of the disease, compared to healthy individuals from the endemic region and a non-endemic area. Methods: We characterized the B-cell repertoire in i) nontreated patients (n=5); ii) patients under immunosuppressive treatment (n=5); iii) patients in remission without treatment (n=6); and two control groups iv) from the endemic (n=6) and v) non-endemic areas in Brazil (n=4). We used total RNA extracted from peripheral blood mononuclear cells and performed a comprehensive characterization of the variable region of immunoglobulin heavy chain (IGH) in IgG and IgM using next-generation sequencing. Results: Compared to individuals from a different area, we observed remarkably lower clonotype diversity in the B-cell immune repertoire of patients and controls from the endemic area (p < 0.02), suggesting that the immune repertoire in the endemic area is under geographically specific and intense environmental pressure. Moreover, we observed longer CDR3 sequences in patients, and we identified differential disease-specific usage of IGHV segments, including increased IGHV3-30 and decreased IGHV3-23 in patients with active disease (p < 0.04). Finally, our robust network analysis discovered clusters of CDR3 sequences uniquely observed in patients with PF. Discussion: Our results indicate that environmental factors, in addition to disease state, impact the characteristics of the repertoire. Our findings can be applied to further investigation of the environmental factors that trigger pemphigus and expand the knowledge for identifying new targeted and more effective therapies.

The Second Highest Prevalence of Celiac Disease Worldwide: Genetic and Metabolic Insights in Southern Brazilian Mennonites.

Celiac disease (CD), despite its high morbidity, is an often-underdiagnosed autoimmune enteropathy. Using a modified version of the Brazilian questionnaire of the 2013 National Health Survey, we interviewed 604 Mennonites of Frisian/Flemish origin that have been isolated for 25 generations. A subgroup of 576 participants were screened for IgA autoantibodies in serum, and 391 participants were screened for HLA-DQ2.5/DQ8 subtypes. CD seroprevalence was 1:29 (3.48%, 95% CI = 2.16-5.27%) and biopsy-confirmed CD was 1:75 (1.32%, 95% CI = 0.57-2.59%), which is superior to the highest reported global prevalence (1:100). Half (10/21) of the patients did not suspect the disease. HLA-DQ2.5/DQ8 increased CD susceptibility (OR = 12.13 [95% CI = 1.56-94.20], p = 0.003). The HLA-DQ2.5 carrier frequency was higher in Mennonites than in Brazilians (p = 7 × 10-6). HLA-DQ8 but not HLA-DQ2.5 carrier frequency differed among settlements (p = 0.007) and was higher than in Belgians, a Mennonite ancestral population (p = 1.8 × 10-6), and higher than in Euro-Brazilians (p = 6.5 × 10-6). The glutathione pathway, which prevents reactive oxygen species-causing bowel damage, was altered within the metabolic profiles of untreated CD patients. Those with lower serological positivity clustered with controls presenting close relatives with CD or rheumatoid arthritis. In conclusion, Mennonites have a high CD prevalence with a strong genetic component and altered glutathione metabolism that calls for urgent action to alleviate the burden of comorbidities due to late diagnosis.

Distribution of a novel CYP2C haplotype in Native American populations.

The CYP2C19 gene, located in the CYP2C cluster, encodes the major drug metabolism enzyme CYP2C19. This gene is highly polymorphic and no-function (CYP2C19*2 and CYP2C19*3), reduced function (CYP2C19*9) and increased function (CYP2C19*17) star alleles (haplotypes) are commonly used to predict CYP2C19 metabolic phenotypes. CYP2C19*17 and the genotype-predicted rapid (RM) and ultrarapid (UM) CYP2C19 metabolic phenotypes are absent or rare in several Native American populations. However, discordance between genotype-predicted and pharmacokinetically determined CYP2C19 phenotypes in Native American cohorts have been reported. Recently, a haplotype defined by rs2860840T and rs11188059G alleles in the CYP2C cluster has been shown to encode increased rate of metabolism of the CYP2C19 substrate escitalopram, to a similar extent as CYP2C19*17. We investigated the distribution of the CYP2C:TG haplotype and explored its potential impact on CYP2C19 metabolic activity in Native American populations. The study cohorts included individuals from the One Thousand Genomes Project AMR superpopulation (1 KG_AMR), the Human Genome Diversity Project (HGDP), and from indigenous populations living in Brazil (Kaingang and Guarani). The frequency range of the CYP2C:TG haplotype in the study cohorts, 0.469 to 0.598, is considerably higher than in all 1 KG superpopulations (range: 0.014-to 0.340). We suggest that the high frequency of the CYP2C:TG haplotype might contribute to the reported discordance between CYP2C19-predicted and pharmacokinetically verified CYP2C19 metabolic phenotypes in Native American cohorts. However, functional studies involving genotypic correlations with pharmacokinetic parameters are warranted to ascertain the importance of the CYP2C:TG haplotype.

Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus-A Possible Epigenetic Effect in the Autoimmune Response.

Pemphigus foliaceus (PF) is an autoimmune skin blistering disease characterized by antidesmoglein-1 IgG production, with an endemic form (EPF) in Brazil. Genetic and epigenetic factors have been associated with EPF, but its etiology is still not fully understood. To evaluate the genetic association of histone (de)acetylation-related genes with EPF susceptibility, we evaluated 785 polymorphisms from 144 genes, for 227 EPF patients and 194 controls. Carriers of HDAC4_rs4852054*A were more susceptible (OR = 1.79, p = 0.0038), whereas those with GSE1_rs13339618*A (OR = 0.57, p = 0.0011) and homozygotes for PHF21A_rs4756055*A (OR = 0.39, p = 0.0006) were less susceptible to EPF. These variants were not associated with sporadic PF (SPF) in German samples of 75 SPF patients and 150 controls, possibly reflecting differences in SPF and EPF pathophysiology. We further evaluated the expression of histone (de)acetylation-related genes in CD4+ T lymphocytes, using RNAseq. In these cells, we found a higher expression of KAT2B, PHF20, and ZEB2 and lower expression of KAT14 and JAD1 in patients with active EPF without treatment compared to controls from endemic regions. The encoded proteins cause epigenetic modifications related to immune cell differentiation and cell death, possibly affecting the immune response in patients with PF.

Cross-sectional study for COVID-19-related mortality predictors in a Brazilian state-wide landscape: the role of demographic factors, symptoms and comorbidities.

OBJECTIVE: The Brazilian state of Paraná has suffered from COVID-19 effects, understanding predictors of increased mortality in health system interventions prevent hospitalisation of patients. We selected the best models to evaluate the association of death with demographic characteristics, symptoms and comorbidities based on three levels of clinical severity for COVID-19: non-hospitalised, hospitalised non-ICU ward and ICU ward. DESIGN: Cross-sectional survey using binomial mixed models. SETTING: COVID-19-positive cases diagnosed by reverse transcription-PCR of municipalities located in Paraná State. PATIENTS: Cases of anonymous datasets of electronic medical records from 1 April 2020 to 31 December 2020. PRIMARY AND SECONDARY OUTCOME MEASURES: The best prediction factors were chosen based on criteria after a stepwise analysis using multicollinearity measure, lower Akaike information criterion and goodness-of-fit χ2 tests from univariate to multivariate contexts. RESULTS: Male sex was associated with increased mortality among non-hospitalised patients (OR 1.76, 95% CI 1.47 to 2.11) and non-ICU patients (OR 1.22, 95% CI 1.05 to 1.43) for symptoms and for comorbidities (OR 1.89, 95% CI 1.59 to 2.25, and OR 1.30, 95% CI 1.11 to 1.52, respectively). Higher mortality occurred in patients older than 35 years in non-hospitalised (for symptoms: OR 4.05, 95% CI 1.55 to 10.54; and for comorbidities: OR 3.00, 95% CI 1.24 to 7.27) and in hospitalised over 40 years (for symptoms: OR 2.72, 95% CI 1.08 to 6.87; and for comorbidities: OR 2.66, 95% CI 1.22 to 5.79). Dyspnoea was associated with increased mortality in non-hospitalised (OR 4.14, 95% CI 3.45 to 4.96), non-ICU (OR 2.41, 95% CI 2.04 to 2.84) and ICU (OR 1.38, 95% CI 1.10 to 1.72) patients. Neurological disorders (OR 2.16, 95% CI 1.35 to 3.46), neoplastic (OR 3.22, 95% CI 1.75 to 5.93) and kidney diseases (OR 2.13, 95% CI 1.36 to 3.35) showed the majority of increased mortality for ICU as well in the three levels of severity jointly with heart disease, diabetes and CPOD. CONCLUSIONS: These findings highlight the importance of the predictor's assessment for the implementation of public healthcare policy in response to the COVID-19 pandemic, mainly to understand how non-pharmaceutical measures could mitigate the virus impact over the population.

Single Nucleotide Variants as Proxies for HLA-A*31:01 in Native American Populations.

Carbamazepine triggers dermatologic hypersensitivity reactions, associated with specific human leukocyte antigens (HLAs), especially HLA-B*15:02 and HLA-A*31:01. Previous efforts to identify single nucleotide variants (SNVs) with high predictive value as HLA proxies, revealed that rs1061235 and rs17179220 fulfill these requirements for HLA-A*31:01 in some but not all populations. Herein we explored the predictive performance of rs1061235 and rs17179220 as HLA-A*31:01 tags in populations of Native American ancestry, which are largely underrepresented in pharmacogenomic studies. The study cohorts comprised the overall Admixed American superpopulation of the 1000 Genomes Project (1 KG_AMR), a subcohort of individuals with predominant Native American ancestry (1 KG_NAT), the Native American population of the Human Genome Diversity Project (HGDP), plus Kaingang (KRC) and Guarani (GRC and GKW) adults from indigenous reservation areas in Brazil. The diversity of cohorts is reflected in the range of frequencies of HLA-A*31:01 (0.02-0.65), rs1061235 (0.03-0.13) and rs17179220 (0.12-0.66), as well as in the predictive performance of these SNVs as HLA-A*31:01 proxies. NPV (negative predictive value), the metric of primary interest for pharmacogenetic-informed carbamazepine prescription was maximal (NPV = 1.0) for both SNVs in 1 KG_AMR and 1 KG_NAT, for rs17179220, but not rs1061235 (NPV = 0.91) in HGDP, and for rs17179220 in GWK, but not GRC (NPV = 0.88) or KRC (NPV = 0.80). Collectively, the data support the notion that rs1061235 and rs17179220 are not optimal proxies for HLA-A*31:01 across populations of Native American ancestry.

Genetic Associations and Differential mRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus.

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

Population structure and forensic genetic analyses in Guarani and Kaingang Amerindian populations from Brazil.

When DNA profile comparisons between a crime scene trace and a reference sample generate correspondence, the match probability has to be estimated, so that evaluation of the strength of the forensic DNA evidence can be made. The random match probability estimations require information on allele frequencies and an adjustment factor, referred to as theta (θ) or Fst, a co-ancestry correction factor for subpopulation effects. The θ value has been standardized for urban and isolated populations, but inconsistencies have been reported when it is specifically calculated for smaller and isolated populations, including Amerindian populations. Notably, attempts to characterize forensic markers of these minor populations have been extensively limited and more conservative estimates of the correction factor may be generated for each of them. Therefore, we estimate allele frequencies of 21 autosomal STR markers used for forensic testing and calculated relevant forensic parameters for the set. In addition, we featured the possible structure of five Brazilian Amerindian populations that have been genetically isolated for centuries so we could obtain the appropriate θ value for them. The sample consisted of 319 individuals: (1) 121 Kaingang, from two communities: Ivaí (KIV=61) and Rio das Cobras (KRC=60); and (2) 198 Guaranis from three communities: Mbya from Rio das Cobras (GRC=51), Guarani Ñandeva (GND=71) and Guarani Kaiowá (GKW=76). Between Guarani populations low (Rst=0.0402, p < 10-4) to high (Rst=0.1557, p < 10-5) differentiation was found. Regarding Guarani and Kaingang populations, intermediate (Rst=0.0590, p < 10-5) to high (Rst=0.1604, p < 10-5) differentiation was found. The two Kaingang populations showed very low differentiation between them (Rst=0.0017, p = 0.27), which justifies the union of both genetic data for forensic databases and calculations. The combined power of discrimination (PD) and the combined power of exclusion (PE) were calculated for each population, demonstrating the usefulness of this set of markers in forensic and kinship analysis regarding these populations. Considering the demographic heterogeneity of Amerindian populations in general, the Fst mean value (0.03) was evaluated regarding 43 different indigenous populations from the Americas, including Guaranis and Kaingangs. This result confirms the adequacy of the standardized θ value for the forensic random match probability estimations involving Amerindian populations.

Uniparental markers reveal new insights on subcontinental ancestry and sex-biased admixture in Brazil.

The Brazilian population is a product of asymmetric admixture among European men and Amerindian and African women. However, Brazilian subcontinental ancestry is scarcely documented, especially regarding its African roots. Here, we aimed to unveil the uniparental continental and subcontinental contributions from distinct Brazilian regions, including South (n = 43), Southeast (n = 71), the poorly genetically characterized Central-Western region (n = 323), and a subset of unique Brazilian Amerindians (n = 24), in the context of their genome-wide ancestral contributions. The overwhelming majority of European Y haplogroups (85%) contrast sharply with the predominant African and Amerindian mtDNA haplogroups (73.2%) in admixed populations, whereas in Amerindians, non-Native haplogroups could only be detected through the paternal line. Our in-depth investigation of uniparental markers showed signals of an Andean and Central-Brazilian Amerindian maternal contribution to Southeastern and Central-Western Brazil (83.1 ± 2.1% and 56.9 ± 0.2%, respectively), the last having the highest paternal Amerindian ancestry yet described for an admixed Brazilian region (9.7%) and contrasting with higher Southern-Brazilian Amerindian contribution to Southern Brazil (59.6 ± 1%). Unlike the higher African Bantu contribution previously reported for the South and Southeast, a relevant Western African non-Bantu contribution was detected in those regions (85.7 ± 5% and 71.8 ± 10.8% respectively). In contrast, a higher Bantu contribution was described for the first time in the Central-West (64.8 ± 1.3% maternal and 86.9 ± 9.6% paternal). We observed sex-biased signatures consistent with the historically recorded Brazilian colonization and added new insights in the subcontinental maternal ancestry of Brazilians from regions never studied at this level.

Extracellular vesicle microRNAs in celiac disease patients under a gluten-free diet, and in lactose intolerant individuals.

BACKGROUND: Celiac disease (CD) is an autoimmune disorder triggered by an abnormal immunological response to gluten ingestion and is associated with deregulated expression of cellular microRNAs (miRNAs) of the gut mucosa. It is frequently misdiagnosed as lactose intolerance (LI) due to symptom resemblance. Microvilli loss may be counteracted by a rigorous gluten-free diet (GFD). AIMS: To identify altered extracellular vesicle miRNAs from plasma among CD patients on GFD (n=34), lactose intolerant individuals on restrictive diet (n=14) and controls (n=23), and to predict biological pathways in which these altered miRNAs may play a part. METHODS: Five different small RNA samples of each group were pooled twice and then screened by new-generation sequencing. Four miRNAs were selected to be quantified by RT-qPCR in the entire sample. RESULTS: The levels of four miRNAs - miR-99b-3p, miR-197-3p, miR-223-3p, and miR-374b-5p - differed between CD patients and controls (P<0.05). Apart from miR-223-3p, all these miRNAs tended to have altered levels also between LI and controls (P<0.10). The results for miR-99b-3p and miR-197-3p between CD and controls were confirmed by RT-qPCR, which also indicated different levels of miR-99b-3p and miR-374b-5p between CD-associated LI and LI (P<0.05). CONCLUSIONS: These miRNAs may have targets that affect cell death, cell communication, adhesion, and inflammation modulation pathways. Hence, altered miRNA levels could be associated with CD-related aspects and gut mucosa recovery.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas.

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

An immunogenetic view of COVID-19.

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

GWAS-Top Polymorphisms Associated With Late-Onset Alzheimer Disease in Brazil: Pointing Out Possible New Culprits Among Non-Coding RNAs.

Several genome-wide association studies (GWAS) have been carried out with late-onset Alzheimer's disease (LOAD), mainly in European and Asian populations. Different polymorphisms were associated, but several of them without a functional explanation. GWAS are fundamental for identifying loci associated with diseases, although they often do not point to causal polymorphisms. In this sense, functional investigations are a fundamental tool for discovering causality, although the failure of this validation does not necessarily indicate a non-causality. Furthermore, the allele frequency of associated genetic variants may vary widely between populations, requiring replication of these associations in other ethnicities. In this sense, our study sought to replicate in 150 AD patients and 114 elderly controls from the South Brazilian population 18 single-nucleotide polymorphisms (SNPs) associated with AD in European GWAS, with further functional investigation using bioinformatic tools for the associated SNPs. Of the 18 SNPs investigated, only four were associated in our population: rs769449 (APOE), rs10838725 (CELF1), rs6733839, and rs744373 (BIN1-CYP27C1). We identified 54 variants in linkage disequilibrium (LD) with the associated SNPs, most of which act as expression or splicing quantitative trait loci (eQTLs/sQTLs) in genes previously associated with AD or with a possible functional role in the disease, such as CELF1, MADD, MYBPC3, NR1H3, NUP160, SPI1, and TOMM40. Interestingly, eight of these variants are located within long non-coding RNA (lncRNA) genes that have not been previously investigated regarding AD. Some of these polymorphisms can result in changes in these lncRNAs' secondary structures, leading to either loss or gain of microRNA (miRNA)-binding sites, deregulating downstream pathways. Our pioneering work not only replicated LOAD association with polymorphisms not yet associated in the Brazilian population but also identified six possible lncRNAs that may interfere in LOAD development. The results lead us to emphasize the importance of functional exploration of associations found in large-scale association studies in different populations to base personalized and inclusive medicine in the future.

Genetic association and differential expression of HLAComplexGroup lncRNAs in pemphigus.

BACKGROUND: Pemphigus is a group of bullous diseases characterized by acantholysis and skin blisters. As for other autoimmune diseases, the strongest genetic associations found so far for pemphigus foliaceus (PF) and vulgaris (PV) are with alleles of HLA genes. However, apart from protein-coding genes, the MHC region includes a set of poorly explored long non-coding RNA (lncRNA) genes, the HLA complex group (HCG). OBJECTIVES: To investigate if HCG lncRNA alleles are associated with pemphigus susceptibility. METHODS AND RESULTS: We analyzed SNPs in 13 HCG lncRNA genes, both in PV (Germany: 241 patients; 1,188 controls) and endemic PF (Brazil: 227 patients; 194 controls), applying multivariate logistic regression. We found 55 associations with PV (pcorr < 0.01) and nine with endemic PF (pcorr < 0.05), the majority located in TSBP1-AS1 (which includes HCG23) and HCG27 lncRNA genes, independently of HLA alleles previously associated with pemphigus. The association of TSBP1-AS1 rs3129949*A allele was further replicated in sporadic PF (p = 0.027, OR = 0.054; 75 patients and 150 controls, all from Germany). Next, we evaluated the expression levels of TSBP1-AS1, TSBP1, HCG23, and HCG27 in blood mononuclear cells of Brazilian patients and controls. HCG27 was upregulated in endemic PF (p = 0.035, log2 FC = 1.3), while TSBP1-AS1 was downregulated in PV (p = 0.029, log2 FC = -1.29). The same expression patterns were also seen in cultured keratinocytes stimulated with IgG antibodies from patients and controls from Germany. TSBP1 mRNA levels were also decreased in endemic PF blood cells (p = 0.042, log2 FC = -2.14). TSBP1-AS1 and HCG27 were also observed downregulated in CD19+ cells of endemic PF (p < 0.01, log2 FC = -0.226 and -0.46 respectively). CONCLUSIONS: HCG lncRNAs are associated with susceptibility to pemphigus, being TSBP1-AS1 and HCG27 also differentially expressed in distinct cell populations. These results suggest a role for HCG lncRNAs in pemphigus autoimmunity.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans.

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Blood groups in Native Americans: a look beyond ABO and Rh.

The study presents comparisons between blood group frequencies beyond ABO and Rh blood systems in Native American populations and previously published data from Brazilian blood donors. The frequencies of Diego (c.2561C>T, rs2285644), Kell (c.578C>T, rs8176058), Duffy (c.125A>G, rs12075, c.1-67T>C, rs2814778) and Kidd (c.838A>G, rs1058396) variants in Kaingang (n=72) and Guarani (n=234) populations from Brazil (1990-2000) were obtained and compared with data from these populations sampled during the 1960s and with individuals of different Brazilian regions. Data showed high frequencies of DI*01 and FY*01 alleles: 11.8% and 57.6% in Kaingang and 6.8% and 75.7% in Guarani groups, respectively. The main results indicated: (1) reduction in genetic distance over time of Kaingang and Guarani in relation to other Brazilian populations is suggestive of ongoing admixture; (2) significant differences in some frequencies of blood group markers (especially Diego, Kidd and Duffy) in relation to Native Americans and individuals from different geographical regions of Brazil. Our study shows that the frequency of red blood cell polymorphisms in two Native American groups is very different from that of blood donors, when we evaluated blood groups different from ABO and Rh systems, suggesting that a better ethnic characterization of blood unit receptors is necessary.

Genetic variability of immune-related lncRNAs: polymorphisms in LINC-PINT and LY86-AS1 are associated with pemphigus foliaceus susceptibility.

Pemphigus foliaceus (PF) is an autoimmune blistering disease of the skin, clinically characterized by erosions and, histopathologically, by acantholysis. PF is endemic in the Brazilian Central-Western region. Numerous single nucleotide polymorphisms (SNPs) have been shown to affect the susceptibility for PF, including SNPs at long non-coding RNA (lncRNA) genes, which are known to participate in many physiological and pathogenic processes, such as autoimmunity. Here, we investigated whether the genetic variation of immune-related lncRNA genes affects the risk for endemic and sporadic forms of PF. We analysed 692 novel SNPs for PF from 135 immune-related lncRNA genes in 227 endemic PF patients and 194 controls. The SNPs were genotyped by Illumina microarray and analysed by applying logistic regression at additive model, with correction for sex and population structure. Six associated SNPs were also evaluated in an independent German cohort of 76 sporadic PF patients and 150 controls. Further, we measured the expression levels of two associated lncRNA genes (LINC-PINT and LY86-AS1) by quantitative PCR, stratified by genotypes, in peripheral blood mononuclear cells of healthy subjects. We found 27 SNPs in 11 lncRNA genes associated with endemic PF (p < .05 without overlapping with protein-coding genes). Among them, the LINC-PINT SNP rs10228040*A (OR = 1.47, p = .012) was also associated with increased susceptibility for sporadic PF (OR = 2.28, p = .002). Moreover, the A+ carriers of LY86-AS1*rs12192707 mark lowest LY86-AS1 RNA levels, which might be associated with a decreasing autoimmune response. Our results suggest a critical role of lncRNA variants in immunopathogenesis of both PF endemic and sporadic forms.

Variation in genes implicated in B-cell development and antibody production affects susceptibility to pemphigus.

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease characterized by the presence of pathogenic autoantibodies against desmoglein 1, a component of intercellular desmosome junctions. PF occurs sporadically across the globe and is endemic in some Brazilian regions. Because PF is a B-cell-mediated disease, we aimed to study the impact of variants within genes encoding molecules involved in the different steps of B-cell development and antibody production on the susceptibility of endemic PF. We analysed 3,336 single nucleotide polymorphisms (SNPs) from 167 candidate genes genotyped with Illumina microarray in a cohort of 227 PF patients and 193 controls. After quality control and exclusion of non-informative and redundant SNPs, 607 variants in 149 genes remained in the logistic regression analysis, in which sex and ancestry were included as covariates. Our results revealed 10 SNPs within or nearby 11 genes that were associated with susceptibility to endemic PF (OR >1.56; p < 0.005): rs6657275*G (TGFB2); rs1818545*A (RAG1/RAG2/IFTAP);rs10781530*A (PAXX), rs10870140*G and rs10781522*A (TRAF2); rs535068*A (TNFRSF1B); rs324011*A (STAT6);rs6432018*C (YWHAQ); rs17149161*C (YWHAG); and rs2070729*C (IRF1). Interestingly, these SNPs have been previously associated with differential gene expression, mostly in peripheral blood, in publicly available databases. For the first time, we show that polymorphisms in genes involved in B-cell development and antibody production confer differential susceptibility to endemic PF, and therefore are candidates for possible functional studies to understand immunoglobulin gene rearrangement and its impact on diseases.